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We study irreversible dimer nucleation on top of terraces during epitaxial growth in one and two 
dimensions, for all values of the step-edge barrier. The problem is solved exactly by transforming 
. it into a first passage problem for a random walker in a higher-dimensional space. The spatial 

distribution of nucleation events is shown to differ markedly from the mean-field estimate except 
in the limit of very weak step-edge barriers. The nucleation rate is computed exactly, including 
numerical prefactors. 
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I. INTRODUCTION 



The understanding of how atomistic processes influence morphology at large scales is of fundamental importance 
for controlled growth of crystalline films via deposition techniques. The irreversible nucleation of immobile dimers, 
giving rise to new terraces, is a key process for the growth of a high symmetry surface. In the preceding paper |Q we 
have shown that for the nucleation on top of existing terraces, the usual Mean Field (MF) theory ^,|) is equivalent 
to considering particles as noninteracting: i. e. not feeling each other even if they are on the same lattice site, so 
that they can meet several times before leaving the terrace. Mean Field Theory (MFT) counts all these fictitious 
nucleation events and therefore leads to an overestimate of the nucleation rate ui, that in most cases is a very poor 
approximation of the correct results. For the spatial distribution of nucleation events we have shown in Ref. [Q that 
a substantial discrepancy between mean-field and exact results is expected, because fictitious nucleations beyond the 
first one always dominate in d — 1 and d = 2. 

In this paper we go beyond mean-field theory and present a series of exact results. We calculate the spatial [-P(n)] 
and temporal [Q(t)] distributions of nucleation events. The quantity Q{t) is the probability that two atoms meet a 
time t after deposition of the second atom and it is formally defined in Eq. (j|) . The evaluation of P and Q allows the 
determination of the total probability W that two atoms meet and this allows the exact computation of the nucleation 
rate u). 

j» The solution of the problem is obtained by mapping the diffusion of two particles on a d-dimensional terrace into 

the motion of a single random walker in d! = 2d dimensions. The statistics of meeting events between the two adatoms 
fNj [ (nucleations) is then obtained as the solution of a suitable first passage problem for the d'-dimensional random walker. 

In d = 1 the problem can be treated analytically in full detail, leading to closed form expressions for all quantities of 
interest. In d = 2 one can easily obtain the results numerically, with arbitrary accuracy. 

Results indicate that the spatial distribution of nucleation sites is very different from the mean-field estimate in 
the limit of strong Ehrlich-Schwoebel barriers, both in d = 1 and d = 2. In the opposite limit of zero or weak 
barriers instead, the difference between the mean-field estimate and the exact result is, for reasonable terrace sizes, 
quite surprisingly small. The temporal distribution of nucleation events decays slowly for short times and later 
exponentially. Finally, the calculation of the nucleation rate to is completed by the rigorous determination of the 
O . nucleation probability W = ~}2 n P{n). This confirms that MFT is safely applicable only for weak barriers in d = 2 
and gives the exact expressions for uj that must be used instead of the MF approximate ones. 

The paper is organized as follows. In Sec. |TJ the model for irreversible nucleation is presented and the fundamental 
quantities needed in the rest of the paper are introduced. The method for the solution of the problem is also outlined. 



Secs.ttll and IV are devoted to the presentation of the exact results obtained in d = 1 and d = 2, respectively. In 
Sec. |V| these results are discussed and interpreted in physically intuitive terms. The conclusions and the perspectives 
of this work can be found in Sec. |v|. 

Some of the most important results have been presented previously in Ref. 0. 
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II. THE PROBLEM AND THE METHOD OF SOLUTION 



In this Section we briefly recall the basic concepts of irreversible dimer nucleation along with some results, obtained 
in the first paper, that will be needed in the following. 

We consider particles deposited onto a crystalline terrace of size L, modeled as a discrete lattice (a square lattice 
in d = 2). The flux of particles is uncorrelated, uniform and of intensity F, so that the average interarrival time is 
Tdep = {FL d )~ 1 . Once on the terrace, an adatom hops at rate (At) -1 = 2dD to a randomly chosen nearest neighbor, 
until it either meets another adatom or leaves the terrace. 

This last process can be hindered by the additional Ehrlich-Schwoebel (ES) barrier H reducing interlayer transport 
to a rate 2dD': the ES length £ ES = (-^ — l)ao measures the strength of the barrier (in the following the lattice 
constant eto is used as unit length). 

The average time spent by a single adatom on the terrace is the residence time and depends on L and £ ES , 

r res = (pL + a£ ES )L/D . (1) 

In the limit ^ ES = 0, r res is equal to Tt r , the average time needed by an adatom to reach the terrace boundary. 
Depending on the value of £ ESl three different regimes may occur: i) Zero or weak barriers (r tr ~ T res <C Tdep)\ ii) 
Strong barriers (r tr <C T res <C T^ ep ); iii) Infinite barriers (r tr <C Td ep <C r res ). 

Particles are deposited according to an exponential distribution of interarrival times: P dap (T) = T de\ ex P( — T / T dep)- 
This implies that all quantities should be computed for a generic interarrival time r and then the results should be 
averaged over P dcp (r). However, we have shown in Ref. jl| that this is equivalent to considering two particles deposited 
simultaneously, one with distribution p^ = 1 / L d and the other with an effective distribution 

P?= T T Pl> ( 2 ) 
7~dep t T res 

where p^ is the normalized solution of the discrete stationary diffusion equation in the presence of a constant flux. 
For infinite barriers (regime iii) p^ 1 — p% = 1/L d . For strong but finite barriers (regime ii) pjf = -^^-p^ = T^Ta- ^ n 
the limit of zero or weak barriers (regime i) p°" = -^F^Pn where p^ has a parabolic shape that vanishes at the edges, 
reflecting the presence of absorbing boundaries. 

Nucleation of dimers takes place when particles are on adjacent lattice sites; here we will assume instead that a 
dimer is formed when two particles are on the same site: this avoids useless mathematical complications without 
modifying the physics of the nucleation process. 

The physical quantities we are interested in are P{n), Q(t) and lu. 

• P(n) is the spatial distribution of nucleation events, computed for two adatoms deposited at the same time with 
normalized distributions p s (the first) and p u (the second). p( N \n) is its normalized version. 

• The distribution Q(t) is the probability that a nucleation event occurs at time t, if the two adatoms have been 
deposited at time zero. Q(t) is not considered within the standard mean-field theory. 

P(n) and Q(t) are derived from the same quantity, the probability R(n,t) that a nucleation event occurs on site n 
at time t: 

P(n) = £)fl(n,t), Q(i) = ^i?(M)- (3) 

t n 

We can also define W, the probability that two atoms meet before leaving the terrace: it is clearly related to P(n) 
and Q(t), because 

w=j2rm = J2 p ^ = I2qw • (4) 

n,t n t 

W is equal to 1 for large or infinite barriers (regimes ii and iii) , but it differs from unity in regime i. The normalized 
spatial distribution is clearly P^ N \n) = P(n)/W. 

• The nucleation rate u> is the total number of nucleation events that occur on the whole terrace per unit time. It 
is related to P(n) or Q(t) via W: 

u = FL d ^ W . (5) 

Tdep ~T T res 

P(n) and Q(t) (and then W) depend on the normalized initial distributions for the two adatoms. Therefore, they 
have the same expressions in regime ii and iii, where p^ is just a constant. From Eq. (|5|) instead, one immediately 
realizes that u) has different expressions in each of the three regimes. 
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We also consider an artificial model where adatoms are independent diffusing particles. They do not stop when 
they meet and each encounter is considered as a (fictitious) nucleation. As shown in Ref. |Q, this model gives exactly 
the same results as mean-field theory. 

The computation of the quantities of interest requires the evaluation of R(n, t). Since we consider irreversible dimcr 
formation, R(n, t) is the probability that two particles diffusing on a d-dimcnsional terrace meet for the first time on 
site n at time t. A method for treating the diffusion of two particles is to take their d+d coordinates as the coordinates 
of a single random walker diffusing on a d! = 2<i-dimensional hypercubic terrace. In this picture a nucleation event 
corresponds to the d'— dimensional walker reaching the d-dimensional hyperplane where the coordinates of the two 
particles are equal. The irreversibility of dimer formation implies that an absorbing boundary condition must be 
imposed on this d— dimensional hyperplane. The probability of dimer formation R(n, t) is then given by the probability 
current orthogonal to the hyperplane. 

More specifically, in d — 1 we pass from two walkers of coordinates n and to to a single walker on a square terrace. 
Nucleation occurs when the walker reaches the diagonal of such a terrace (n = m) . In d = 2 we must consider a single 
walker in a four-dimensional space whose coordinates are (n\ 1 mi, n%, m^) and the hyperplane is now a bidimensional 
plane defined by the conditions n\ = 712 and mi = m,2. 

In this way we have reduced the dimer nucleation problem to a first passage problem f6f|. The solution of such a 
problem M is possible analytically in d — 1 (Sec. 



Ill) and numerically in higher dimension (Sec. IV). 



III. RESULTS IN ONE DIMENSION 



When the system is one-dimensional, the two adatoms are mapped into a two-dimensional walker hopping inside a 
square lattice of size L with probability Pm,n{t) to be in site (m, n) at time t. Assuming that, when two adatoms are 
present, one of them, randomly chosen, moves once every time unit, the discrete evolution equation for p m ,n(t) is 

Pm.n(t + 1) = -r[Pm+l,n(*) + Pm-l,n{t) +pm,n+l(*) + Pm,n-l(t)] (6) 

where the discrete time unit corresponds now to a physical time At = l/(2d'D) = l/(AdD). 

The indices m, n vary between 1 and L, but in order to use Eq. ([j|) for all terrace sites it is useful to introduce 
fictitious sites in m = 0, L + 1 and n = 0, L + 1. In this way, boundary conditions are easily written for generic values 
of the ES barrier: p 0>n = api >n , pl+i,u = apL,n, Pmfi = ap m ,i, Pm,L+i = a,Pm,L, where a = 4s/ (1 + 4s)- They apply 
at any time and for all edge sites. There is also an additional boundary condition along the square diagonal 

Pn,n(t) = n= 1, ...,L , (7) 

because the two adatoms stop diffusing when they meet. The initial condition is p m .n(0) = PmPn, but it is also correct 
to write p m , n (0) = p^Pn- ^ n order to obtain a spatial distribution P(n) that is properly symmetrical with respect to 
the center of the terrace we use the symmetrized expression: 

Mo) = |bK+«]- (8) 

The basic quantity we want to compute, the nucleation probability on site n at time t + 1, is 

R(n,t+ 1) = |[p n +l,n(*) +Pn-l,n(t) +Pn,n+l(t) + Pn,n-l(t)] ■ (9) 

In the case of noninteracting particles, the boundary condition along the diagonal is dropped and Eq. @ is replaced 
by 

R(n,t) =p„, n (t) . (10) 

An explicit analytic solution of the problem, both for interacting and noninteracting particles, is possible in the 
limits of zero and infinite ES barriers and will be presented in detail below. As remarked in Sec. |l| for P, Q and W 
only two distinct regimes exist, and 4s = 0,oo are their representative limits. 
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A. Zero barriers (regime i) 



When no ES barrier is present, £ ES — and a = 0. Hence the boundary conditions are simply po,n = PL+i,n — 
Prn,o = Pm,L+i = 0; indicating that edges are absorbing boundaries. In the limit £ ES — the normalized stationary 
distribution is 

P »= L(L + l)(L + 2) ^ L + 1 - n) - 



1. Noninteracting adatoms 



By separating space and time variables in a way perfectly analogous to the treatment of a single particle 0] , we 
find the general solution of Eq. (|^) 



L I 
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kit 
L + l 



L + l 



where the coefficients Bkj are 
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Given the explicit form of p m , n (0), the coefficients Bkj are (see Ref. jjj) 



B kj = Ay A 
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and this allows the evaluation of all the quantities of interest. 
• Nucleation sites. The spatial distribution of nucleation sites is: 



(12) 



(13) 

(14) 
(15) 



Pm (n) = R{n, t) = ^ Pn , n (t) 



(16) 



Its normalized version P„f \n) is plotted in Fig. 0. As proven in Ref. Q within a continuum formalism, it is equal to 
the mean field distribution. This result can be easily proven analytically in a discrete lattice as well. 
• Nucleation times. The distribution of nucleation times is given by: 



Qm(t) = y^Pn.nft) 
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(17) 
(18) 
(19) 



and it is plotted in Fig. ^. To find analytically the behavior of Qmit) for large L, we rewrite it in the following form: 
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L + l 
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kn 
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(20) 
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FIG. 1. Normalized spatial distribution p( N \n) for d = 1 and £ ES = 0. Empty circles are for interacting particles, full circles 
for noninteracting particles (Mean Field Theory). L — 100 in the main part, L = 20 in the inset. 




FIG. 2. The temporal distribution Q(t) for d = 1 and L — 100. From top to bottom, data are for noninteracting particles 
(£es = oo and £ ES = 0) and interacting particles (£ E s = oo and £ E s = 0) . The main part of the figure highlights the power-law 
decay for short times (log-log plot): the solid line goes as t^ 1 ^ 2 . The inset highlights the exponential decay for long times 
(lin-log plot). 
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The coefficients B k k diverge for small k as k 4 (-Bfefe — 192/[7r 4 fc 4 (i+l)(i+2)L 2 ]), so that the dominant contribution 
to the sum for large L comes from the first mode k = 1. Expanding the small argument of the cosine this gives 



VI (t) ~ cxp 



t 



2 V L + 1 



(21) 



which is exactly the exponential decay appearing in Fig. @. 

• Nucleation rate. Let us first compute W NI — P(n) = ^2 t Q{t) that for noninteracting particles is not a probability 
but the total number of times the two adatoms meet before leaving the terrace (W NI can be larger than 1): 



L+l 



B 



kk 



Using the explicit form of Bkk , we obtain 



fc=l 1 - cos (^j) 

- 2 (^)sin 2 (^ 



L 2 (L + l) 2 (L + 2) 



E 



For large L the dominant contribution is provided by the term with k = 1, which gives 



VF NI 



L 2 (L + l) 2 (L + 2) 
Hence the total nucleation rate is 



2(L + 1) 



7T J L 2 (L + 2) 



0.2L . 



w N i ~ 

1~dep 



(L + l) 4 
7T J L 2 (L + 2) 



Using the explicit expression H r res — r fr = L 2 /(12D), and considering only the leading order in L we find 



1 /2\ 6 F 2 L 5 



4 V 7T 



(22) 



(23) 



(24) 



(25) 



(26) 



2. Interacting adatoms 



For interacting adatoms it is possible to take advantage of the noninteracting solution (12) by using a trick: we 



pass from the initial condition p m rl (0) [given in Eq. (jg)] to an auxiliary antisymmetric initial condition 

!Pm,n(0) for m < n 
for m = n (27) 
-Pn,m(0) for TO > 71 

which satisfies the boundary condition p n ^ n — along the diagonal. 

Let us observe that the dynamics given in Eq. (J^j conserves the parity of the spatial distribution, because: 

Pm,n(0) = -pn,m(ty B k j = ~B jk =*> P m ,„{t) = -p„,m{t) ■ (28) 

This means that if we start with an antisymmetric distribution p m ,n(0) = —p n ,m{0), the boundary condition 
Pn,n(t) = is obeyed for all times: the two triangles (to > n) and (m < n) are dynamically disconnected. The 
solution of Eq. (||) is therefore still given by Eq. ( p"2| ) , since the boundary condition is fully taken into account by the 
value of the coefficients Bkj, which depends on the antisymmetric form ofp m>n (0). 

The coefficients Bkj are given by 



m,n— 1 



TO/C7r \ . / 7JJ7T 



G 



In this expression we can decompose the summation Ylm n as E m <n + Em> n i m the latter interchange the dumb 
indices n,m and exploit the antisymmetry of p m ,n(0). We finally obtain: 



B k j - 
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mkn 
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where 



L 



mkir 
L + l 



njn 
L + l 



(31) 



The evaluation of B^j is here less straightforward than in the nonintcracting case. In particular, some sums are 
not easily performed explicitly. This makes difficult the presentation of explicit results. Therefore in the following we 
will present only the general results, leaving the coefficients Bkj indicated. 

• Nucleation sites. Since the regions m < n and m > n are equivalent, the probability of a nucleation event on site 
n at time t + 1 > is given by R(n, t + 1) = ^[p n ^ n+ i(t) + |> n _i jn (t)], while for t = nucleations occur because both 
adatoms are deposited on the same site, i. e. with probability p„, n (0). The spatial distribution of nucleation sites is 
therefore 



_^ r.x. 
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(32) 
(33) 
(34) 



and the normalized distribution pW (n) is presented in Fig. ^. The plot clearly shows that the spatial distribution is 

very similar to the mean-field result P^f \n), although a small discrepancy exists (see the discussion in Sec. VA). 
• Nucleation times. The distribution of nucleation times for interacting particles is 



Q(o) = X>«>«(°) . 

71 = 1 
L 1 

Q(t + 1 > 0) = «[Pri,n+l(<) +Pn-l,n(*)] • 



(35) 
(36) 



Once summed, the two terms in Eq. fl36| ) are equal. Hence, 
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Q(t) is shown in Fig. 0. For short times it decays slowly (as t^ 1 / 2 ), while it goes down exponentially for large 
times. In Appendix [a| we show in detail that the behavior of Q(t) for short and long times can be derived explicitly 
in the case of two adatoms with uniform initial distributions, for which the coefficients Bkj are explicitly known: for 
t < 2L 2 /tt 2 we find 



Q(t) 



Lir 5 / 2 y/t/2 



and for t > 2L 2 /tt 2 



Q(t) 



80 
9vr 2 L 2 



exp 
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(42) 



(43) 



No qualitative change is expected if one atom is initially distributed according to p^ rather than p^: only prefactors 
are expected to be different and this is confirmed by the behavior shown in Fig. ||. 
• Nucleation rate. The probability W of a nucleation event is 



W 



1 L 

7+Y. 



kj ^kj 



(44) 



In Appendix ^ we prove that for large L, W goes to a constant. This constant is found numerically to be roughly 
equal to 0.47. Hence, for large L, the nucleation rate is 



PL 2 F L 

FL- — FLW ~ 0.04 

D D 



(45) 



B. Strong and infinite barriers (regimes ii and lit) 

With infinite ES barriers, £ ES — oo and a = 1. Step edges are perfectly reflecting barriers and boundary conditions 
are: pi. n = po.n, Ph,n = PL+i.n, Pm,i = Pm,o, Pm,L = Pm.L+i- The normalized stationary distribution is simply 
Pn=Pn — l/L, because the distribution of the first adatom is still flat when the second arrives. 

The general solution for a two-dimensional walker is now 
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where 
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X k (n) = tan ( — ) sin ( —j— ) + cos ( —j— 
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(47) 



and the coefficients are 



1 



(48) 



with (Sko is the Kronecker symbol): 



1 + tan^ 



kn 
2L 



(l + S k0 ) 



(49) 



1. Noninteracting adatoms 

The case with noninteracting adatoms is completely trivial for infinite barriers. At any time, p m .n(t) = 1/L 2 so 
that the spatial and temporal distributions of nucleation events are constant. The total number of nucleation events 
W m is clearly infinite. 
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FIG. 3. Normalized spatial distribution P^ N \n) for d = 1 and ^ ES = oo. Empty circles are for interacting particles, full 
circles for noninteracting particles (mean-field approximation). L = 100 in the main part, L — 20 in the inset. 

2. Interacting adatoms 

In a way analogous to the case with zero barriers, we consider antisymmetric initial conditions and we obtain 



B kj = j^- P--n(O) {X k (m)X 3 (n) - (k O j)} = [B< - B<] 



where 



B% = E Pm,n(0)^ fc (m)J0(T 



More explicitly (p m , n (0) = 1/L 2 ) 



n-l 



m= 1 



R ^ n=l 

1 L 

AT, E 



L 2 N k N 3 — 

n-l r 

E 



ITT \ . / TLJ7T 

tan — sin 

2iy i 



. kn\ . ( mkir 

tan — sm 

2LJ \ L 



mkir 



• Nucleation sites. The distribution of nucleation sites is given by 

_^ r.x. 

P(n) = p„,n(0) + - }^\p n<n+ l(t) +p n -l,n(t)} 
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and it is plotted in Fig. |^ [in this case p( N \n) and P(n) coincide, since W = 1]. The distribution has a rounded peak 
in the middle of the terrace and vanishes towards the boundaries. 

The above expression for P(n) is exact, but it is not easy to use in applications. A simpler, approximate, expression is 
therefore highly desirable. In Sec. |V A we show that P(n) is well fitted by a hyperbolic cosine. Up to the normalization 
factor: 



P(n) — cosh(7r) — cosh 



2n 
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(57) 
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FIG. 4. Comparison of the normalized spatial distribution P^ N \n) for d — 1 and / ES = oo (circles) with the approximate 
formula ( |57| ) (solid line). L — 128 (main), L — 16 (inset). 



In Fig. |^ we compare exact and approximate distributions: the agreement is fair already for relatively small sizes 
and very good for large sizes. 

• Nucleation times. As in the case with no barriers, we have 



Q(0) = £>„,„(0) , 

n=l 
L 

Q(t+l>0) = 5^p n _i,n(t) 
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Using the expression for p m ,n(t) 



Q(o) 



i 

L ' 

L-l 



Q(i+1>0) = ^ ByCjy 
fc,j=0 



where the coefficients Cjy are now 
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(62) 



The form of Q(t) is shown in Fig. ^. The decay is the same as for zero barriers: for short times it decays as t~ 1//2 and 
for large times exponentially. Physically intuitive interpretations of these behaviors are discussed in Sec. VB. 
• Nucleation rate. Since W = 1, the nucleation rate in regime ii is 
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while in regime Hi it is simply 



u(L) = FL= 



(64) 
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FIG. 5. Normalized spatial distribution P^ N \n) fJr d = 1 and L — 50. EmrfBy circles are for interacting particles, full circles 
for noninteracting particles (MFT). £ ES = (top left), £ ES = 10 (top right), £ ES = 50 (bottom left), £ ES = 250 (bottom right). 



C. Intermediate barriers 



For intermediate values of the barriers, i. e. values of a between and 1, an explicit analytic solution of the problem 
is not possible, even for noninteracting adatoms. This is a direct consequence of the lack of an explicit solution for 
intermediate barriers even in the case of a single particle (Ref. ji}). Nevertheless the problem can easily be solved 
numerically for any ^ ES , through direct calculation of the dynamical evolution of p min (t), which determines R(n,t) 
and all the quantities of interest. 

The systematic error in the results, due to the integration of Eq. (^]) up to a finite time, is fully negligible for 
realistic values of L: the probability Q(t) that a nucleation occurs at time t decays exponentially for large t and 
consequently the systematic error can easily be made exceedingly small. All numerical results presented in this paper 
can be considered virtually exact. 

As expected, the results for intermediate barriers smoothly interpolate between the two limits of zero or infinite 
barriers, being £ ES / L the only relevant parameter. 

The spatial distribution of nucleation events P^ N \n) is presented in Fig. [j] for L = 50 and several values of £ ES 
(£ ES = 0,10,50,250). Even a small value £ ES /L — 1/5 changes in a notable way the distribution P( N '(n), The 
temporal distribution Q(t) smoothly interpolates between the two limit behaviors presented in Fig. |^. 



IV. RESULTS IN TWO DIMENSIONS 



When the terrace is two-dimensional the motion of two adatoms can be mapped into a four-dimensional problem 
for a single random walker: Pm 1 ,n 1 .m.2,n 2 (t) is the probability of finding one atom on site (mi,ni) and the other in 
(ma, TI2) at time t. Such a probability obeys the equation of motion 

Pmi ,m ,m 2 ,n 2 (t T 1) = - [ Pmi+l,ni,m 2 ,n 2 V") ' Proi-l,ni ,m 2 ,n 2 (t) + (65) 

Pmi ,n\ ,7712 + 1,^2 \r) Pin\ ,n\ ,ni2 — X,ri2 (^) 
Pmi ,m ,ni2 ,^2 + 1 ~"~ Prn 1 , n 1 , rn 2,^2 — 1 (t)} 

with the boundary condition p A+ g = apa, where h is any edge site of the four-dimensional hypercube and (h + S) is 
a nearest neighbour outside the cube. The initial condition is 

Pmi,ni,m2,n 2 (0) = 2 \Pmi,niPm2,ri2 Pm 1 ,niPm2,n 2 \ ' (66) 

where - as usual - n = 1 / L 2 is the uniform initial distribution in two dimensions and p^ n is the normalized 
stationary solution of the discrete diffusion equation in d — 2. 
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FIG. 6. Normalized spatial distribution P^ N \n,n) along the diagonal for d — 2, L = 32 and £ ES = 0. Empty circles are for 
interacting particles, full circles for noninteracting particles. 

A. Zero barriers (regime i) 

1. Noninteracting adatoms 

For noninteracting adatoms the computation of the quantities of interest proceeds along the same lines as in the 
one-dimensional case. The general solution of the equation of motion for the four-dimensional random walker is 



Pm 1 ,ni,m 2 ,n 2 (t) — ^ ] L>k 1 j 1 k 2 j 2 



1 



COS 



L + l 



Jin 
L+l 



k 2 n 
L + l 



L+l 



. ( m 1 k 1 n\ . / niji7r\ . / m 2 fc 2 7r \ . / n 2 j 2 ir 
sm I — — ] sin I — ] sin | — — | sm ' 



\ L + l \L + 1 



L + l 



L + l 



where the coefficients B 



B 



L + l 



mi ,ni ,m2 ,ri2 — 1 



, m 1 k 1 ir\ . / n 1 j 1 n\ . ( m 2 k 2 Tr\ . ( n 2 j 2 ir 
sm [ — — | sm | — | sm | — — | sm 



kijik 2 j2 



L + l J ~"~ V L + 1 J ~~ V L + 1 J \L + 1 
Given the initial condition (|6^), the coefficients Bk 1 j 1 k 2 j 2 are °f the form 

o _ aU aS 

JD k 1 j 1 k 2 j 2 — -"-kiji-n-kaja ' 

where -A^ are the coefficients of the expansion of Pm'n (see Ref. jlj). 

The probability R(m,n,t) of a (fictitious) nucleation event at time t on site (m,n) is given by p m ,n,m,n(t). 
• Nucleation sites. The spatial distribution of nucleation sites is 



(67) 

(68) 
(69) 

(70) 
(71) 

(72) 



P NI (m,n) = R(m, n, t) = ^p m , w (t) 



(73) 



and its normalized version is reported in Fig. |^ [we plot it along the diagonal of the square terrace, P^\n, n)]. 
• Nucleation times. The distribution of nucleation times is 
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FIG. 7. The temporal distribution for d — 2 and L = 40. In the main part 1/Q(t) is plotted vs t to highlight the logarithmic 
decay for short times. From top to bottom data are for noninteracting particles (£ ES = oo and ^ B s = 0) and interacting particles 
(£ E s = oo and ^ ES = 0). The solid line goes as ]n(t/to) [i. e. Q(t) ~ 1/ \n(t/to)]- The inset shows the exponential decay of Q(t) 
for long times. 



l(^) — ^ P7n,n,m,n(ty 
m.n—1 
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ki-,ji-k 2 ,j2 = l 

kiir 
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mkiir 
L+l 



cos 
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sm 



L + l 
L+l 
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njiir 



J2TT 

L + l 
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2_j Bktjtktjt S _ 



kl,jl=l 



L + l 

klTT 

L + l 



sm 



L + l 



L + l 



(74) 
(75) 

(76) 
(77) 

(78) 



and it is plotted in Fig. 0. For large times, only the most slowly decaying mode fci = = 1 contributes to the sum, 
yielding 



Qm(t) ~ exp 



t In cos 



L + l 



exp 



t f 7T 



2 V L + 1 



Nucleation rate. The total number of times the two adatoms meet before leaving the terrace is 



W m = Pm{rn,n) 



L + l 



2 L 



E 



B 



m,n—l 

For large values of L only the mode (1,1,1,1) dominates the sum, so 



cos(^) + cos 



(79) 



(80) 



L + l 



C0S (lTi) + cos (lTi) 



_2_ ( L + l 

tt 2 { L 



and the total nucleation rate is 



L + l 

L 



(81) 



(82) 
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FIG. 8. Normalized spatial distribution p( N \n,ri) along the diagonal for d — 2, L — 32 and £ ES = oo. Empty circles are for 
interacting particles, full circles for noninteracting particles. 



Using the expression r res ~ ^ derived in Ref. |IJ] , the leading term in L is 



64 F 2 L 6 
^ D 



(83) 



Interacting adatoms 

At odds with what happens for the one-dimensional case, the trick of using an initial condition antisymmetric with 
respect to particle interchange can not be used in two dimensions for taking into account the interaction between 
particles. The physical reason is that in two dimensions two particles can swap their position without meeting. 
As a consequence, the configuration space can not be split into two dynamically disconnected regions, because the 
condition p m ,n,m,n(t) = holds on a two dimensional plane that does not divide the four-dimensional configuration 
space into separate domains. Hence it is not possible to implement the additional boundary condition p m ,n,m,n{t) = 
by choosing the initial condition to be antisymmetric. We have not been able to overcome this problem analytically 
and therefore for the interacting case we res ort to the numerical solution of Eq. (|65|) , which is easily performed, and 
gives virtually exact results (see Sec. HIC). The results for p( N \n,n) and Q(t) are presented in Figs. [| and [t], 



respectively. The spatial distribution as given by MF theory agrees with exact results even better than in d = 1; the 
short time decay for Q(t) does not follow a power law but rather a logarithmic one [Q(t) ~ 1/ \n(t/to)}. 
For what concerns the total nucleation rate we find numerically W ~ 0.25/ ln(L/1.3) and this implies 

^ - 00 4mWL3)- ^ 



B. Strong and infinite barriers (regimes n and in) 

The results for Q(t) are presented in Fig. |^ and those for P( N )(m, n) in Figs. ||,||. The spatial distribution along the 
diagonal (Fig. @) behaves much in the same way as in d = 1; a qualitatively similar behavior is found along different 
directions (Fig^g). A deeper analysis is deferred to Sec. |V A\ 

The total number of nucleation events W is clearly 1 and this implies, in regime ii 

t F 2 T, 5 P 

u = FL 2 -^= ES (85) 

Tdep 4£> 

and in regime Hi 
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FIG. 9. Normalized spatial distribution P'"' for d — 2, L = 2>2 and £es = 0. Full circles are for noninteracting particles. 
Empty symbols are for interacting particles: along the diagonal (circles), along one edge (squares), and in the middle, parallel 
to one edge (diamonds). 
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FIG. 10. Normalized spatial distribution P'"' (n, X) along the diagonal for d Q 2 and L = 20. Empty circles are for interacting 
particles, full circles for noninteracting particles. £ ES = (top left), £ ES — 4 (top right), £ ES = 20 (bottom left), ^ ES = 100 
(bottom right). 
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7~dep 



(86) 



For intermediate barriers, the results for P^ N \n,n) are presented in Fig. [T(]. As in d = 1, the spatial distribution 
interpolates between the two limits of zero and infinite barriers. As already remarked for the one-dimensional case, a 
relatively small ES barrier (£ ES /L = 1/5) affects quite dramatically the spatial distribution. 



V. DISCUSSION OF THE RESULTS 
A. The spatial distribution 



The form of the spatial distribution of nucleation sites has been presented in the preceding two Sections both in 
one and two dimensions and for all values of the ES barrier. Some remarks are in order. 

As expected pL we find that the mean- field assumption for the distribution of nucleation sites is in general not 
exact both in one and in two dimensions, for all values of £ ES and for all L. The origin of the discrepancy between 
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the exact form of P(n) and the MF counterpart is clear: the mean-field approximation is equivalent to considering 
particles as noninteracting, i. e., taking into account not only the first nucleation event between the particles, but 
also all subsequent encounters between them that would occur should they keep diffusing after meeting. 

Although not exactly the same, the mean-field distribution is however a very good approximation of the true spatial 
distribution, for zero or weak ES barriers, particularly in d = 2 but also in d = 1. This result is somewhat striking, if 
we consider that the ratio W NI /W is proportional to L in d = 1 [Eq. (§4))] and to InL in d = 2 [Eq. (pl|)1. Hence the 
relative weight of successive nucleations diverges for growing L; nevertheless P N i(n) is very close to P(n), indicating 
that the distribution of all nucleation events following the first one is very similar to the distribution of the first. 

Things are radically different for large ES barriers. In this case the discrepancy between MF and true distributions 
is remarkable. Also this result is somewhat counterintuitive. Particles are distributed uniformly at the beginning and 
each of them would remain like that forever in the absence of the other: this is the reason why the spatial distribution 
P NI for noninteracting particles is uniform. The interaction between particles breaks this uniformity. Consider for 
example the one-dimensional case. The nucleation probability on site n close to the center of the terrace, is the sum 
of the statistical weight of all pairs of random walks (one for each particle) with the constraint that they intersect for 
the first time in n. If the site n is close to an edge, one of these walks is reflected by the boundary and the weight of 
walks intersecting for the first time in n is strongly reduced. 

An 'entropic' mechanism is present for weak barriers as well: in this case nucleation close to an edge is made 
difficult by adsorbing boundaries, which reduce the probability to find an atom close to the steps. For weak barriers 
P(n) is peaked around the middle of the terrace also because the initial distribution for one atom is not uniform, but 
parabolic. 

Notice however that for infinite barriers the mean-field distribution, which includes the contribution of successive 
encounters, is completely flat. This indicates that, even if it is relatively unusual for particles to meet close to edges, 
once this happens they tend to meet there several times and this restores uniformity in the distribution of all nucleation 
events. 

In Eq. ( |57| ) we have proposed an approximate expression for the distribution P(n) in the limit of infinite ES barriers. 
It has been derived assuming a behavior as an hyperbolic cosine 

P(n) = ao[ai — cosh(a 2 n — 03)] (87) 

and imposing that P(n) is symmetrical with respect the center of the terrace (03/0,2 = (L + l)/2) and that P(0) = 
P(L + 1) = [ai = cosh(a 3 )]. The former condition is obvious, and the latter derives from the numerical evidence 
that P(l)/P (£±1) goes to zero for increasing L. Once both conditions have been imposed we obtain: 



P(n) = a 



cosho.-j cosh [ — ^ — 03 



(88) 



There is only one fit parameter, a 3 , because ao is constrained by the normalization condition for P(n). From a 
nonlinear curve fitting for relatively small L we can extrapolate that a 3 (L) tends to a constant value of order 3.11 as 
L grows. We have, somewhat arbitrarily, set 03(00) equal to 71". 



B. The temporal distribution 



The results for the temporal distribution of nucleation events Q(t) show in all cases a slow decay for short times 
(as a power- law in d = 1, logarithmic in d = 2) followed by an exponential decrease for larger times. This behavior 
has been obtained by solving exactly analytically or numerically, the evolution equation for the particles on the 
terrace. Its physical meaning is clarified further by rederiving these results by means of more transparent but less 
rigorous arguments: the decay for short times is interpreted in terms of first passage properties of random walks in 
an unbounded space; the long time decay is the combined effect of the exponentially decreasing probability that both 
particles are still on the terrace at time t and the probability that they have not yet met. 

Let us first discuss the behavior at short times and consider the relative coordinate of the two particles as the 
coordinate of a fictitious particle C: nucleation occurs when C reaches the origin. The initial spatial distribution 
probability for C [pcM] is a function of pa and ps complicated by the presence of boundaries. However, we are 
interested in the behavior for short times, i. e. times such that particles are not affected by the presence of terrace 
edges. Therefore we can assume an initial spatial distribution p c {r) uniform in a region of linear size L around the 
origin (p c = l/L d ) and zero outside. The irrelevance of boundaries in the short time regime is confirmed by Figs. ^ 
and0: Q(t) has the same behavior, independently of step-edge barriers. 
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We now define F(r, t) as the first passage probability in r at time t starting from the origin at time zero. The 
probability that atom A, leaving from r at t = arrives for the first time in the origin at time t is clearly F{—r, t), 
so that 

Q(*)=$>o(r)F(-r,t) • (89) 

r 

Let us also define P(r, t) as the probability that a particle is in r at time t, being at the origin at time zero. At t = 
we have P(r, 0) = S rt o and F(r, 0) = 0. P(r, t) and F(r, t) are connected by JsJ 



P(r,f) = ^F(r,T)P(0,f -r) . (90) 

T=0 

We write Eq. ( |90| ) for spatial argument — r, multiply both sides by pa(f)' 

t 

p c (r)P(-r,t)=J2po(r)F(-r,T)P(0 ) t-T) (91) 

T=0 

and sum over r: 

t 

£ Pc (r)P(-r, t)=J2 Q(r)P(0, t - r) . (92) 

r r=0 

At short times P(r, t) is negligible in the region where p c {r) vanishes. Therefore we can take p c out of the summation 
and use the normalization of P(— r, i), obtaining: 

1 * 

_ = ^Q(r)P(0,t-r) . (93) 

T=0 

In d = 1, we pass to the continuuum in time [P(0, t — t) = l/(i — t) 1 / 2 ], 

(94, 



i A (f-T)'/2 

and setting t = ts, we obtain 



T = / *>7^i7S . (95) 



i 7o (l-.) 1 / 2 



which implies Q(t) ~ f 1 / 2 . 

In two dimensions we separate the term r = t in Eq. ( |93"|) 



1 - f* frQQ+Qfl ( 96) 



L 2 In t-T 



n 



1 ^ ds^- + Q(t) (97) 
1 — s 



' 1 ~ 1/t ds 

Q(t) I t^- + Q(t) (98) 



1 - s 

~Q(i)[l + lni] (99) 

and we obtain Q(i) ~ 1/(1 + hit). 

In conclusion, at short times Q(t) decays as a power law [Q(t) ~ l/Vf] in ei = 1 and logarithmically [Q(t) ~ 
l/ln(t/t )] in d = 2. 

Let us consider now the behavior for long times. The probability that a single adatom remains on the terrace up 
to time t is (see Ref. [jjj) 5(f) ~ exp(— agt). In <i = 1 one has as = \ ( jxj J for zero barriers and as = for infinite 
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TABLE I. Value of the nucleation rate w, including the correct prefactors. In the regime Hi of infinite barriers, lo = (FL d ) 



barriers. It is better to introduce a continuous time notation. The time step for a single particle is At = l/(2dD), 
while for two particles diffusing on the same terrace is At — l/(4dD), so that 

S(t) ~ exp(-2asdDt). (100) 

The probability that two adatoms meet at time t decays as 

Q(t) ~ exp(-4a Q dDt). (101) 



In Sec. Ill we determined that for noninteracting particles ag = | (^XTrJ f° r zero barriers and ctQ = for infinite 

barriers, while for interacting particles ojq = | ( j^pjj for zero barriers and oiq — j (j;) 2 for infinite barriers. 

All these findings are simply rationalized by the following argument. We define G(t) as the probability that two 
adatoms confined on the terrace meet for the first time at time t: it is therefore equal to Q(t) in the limit of infinite 
barriers. For long times, 

G(t) ~ exp(-4a G dDt) . (102) 

We claim that, for interacting adatoms, Q(t) is given by the probability that each of the two adatoms is still on the 
terrace times the probability that they meet for the first time at time t: 

Q(t) ~ S 2 (t)G(t) =► a Q =a s + a G . (103) 

In the noninteracting case, clearly G(t) does not play any role. Then 

Q(t) - S 2 {t) => a Q =a s . (104) 



If we neglect the differences between L and L + 1 at the denominators of «q, the relations ( |103| ) and (104) are both 
verified in the limits £ ES = and I 

In d = 2 the value of ag is not known analytically. However relations (103 104) have been verified numerically. 



C. The nucleation rate 



In this paper we have computed exactly the scaling of the nucleation rate in all regimes, in d = 1 and d = 2, both 
for the noninteracting (mean-field) and the interacting case. The results confirm those of Ref. [Q, where the rigorous 
calculation of W was lacking. Mean field theory overestimates the nucleation rate by a factor that scales, in regime 
i (zero or weak barriers), as L in d = 1 and hiL in d = 2. In the limit of strong barriers (regime ii) the error scales 
as £ ES in d = 1 and as l ES /L in d = 2. Notice that the latter is a large quantity, since in this regime £ ES ^S> L. For 
infinite barriers (regime Hi) the mean-field picture trivially breaks down. Hence mean-field theory is generally strongly 
inaccurate, except in two dimensions for weak barriers; however, even in this case logarithmic corrections render lo mf 
not completely reliable. 

Our treatment allows the evaluation not only of exponents, but also of prefactors. In particular, this is performed 
analytically for zero or strong barriers in d = 1 [Eqs. ( [l5| ) and (|6^)] and for strong barriers in d = 2 [Eq. (p5|)], while 
for d = 2 and no barriers we have evaluated the prefactor numerically. For reference we report in Table | the value of 
the nucleation rate in the different cases. 

Finally, we want to remark that in d — 1 for zero barriers, not only the asymptotic behavior for large L, but the 
exact value of lu for any L, can be determined analytically. One just needs to perform the sum of L 2 terms [Eq. (0)]. 
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VI. CONCLUSIONS 



In this paper and in the preceding one we have presented a rather complete study of the problem of irreversible 
dimer nucleation on top of terraces during epitaxial growth. We have analyzed in detail the mean-field approach to 
this problem, identified its weaknesses and provided a physical interpretation for them. Then we have solved the 
problem, by analytical means or (when needed) numerically In this way we have derived exact results for the spatial 
and temporal distributions of nucleation events and for the total nucleation rate. 

We believe that these results provide a relevant contribution to the investigation of crystal growth both from the 
experimental and the theoretical point of view. 

The dependence of the nucleation rate uj on the terrace size L and the ES length £ ES is a crucial piece of information 
for the interpretation of experimental results, for example the evaluation of the Ehrlich-Schwocbel barrier. The mean- 
field approximation has been widely used so far: as already pointed out ||, this introduces a systematic underestimate 
of the strength of the ES barrier. The exact expressions for lo, derived in this work, must replace the MF approximate 
formulas for a correct interpretation of experimental data. 

From the theoretical point of view, also the spatial distribution plays an important role. Sometimes, 'mesoscopic' 
models are used to describe the growth process in the submonolayer regime [[l0| or in the multilayer regime |TT| . 
The rule for dimer formation must be supplemented with the spatial distribution P(n) of nucleation sites: as we 
have argued, if additional step-edge barriers are not negligible, exact results are completely different from mean-field 
predictions. 

Let us finally mention some possible extensions of the present work. In this paper and in the previous one, we have 
discussed irreversible nucleation on top of compact terraces: it is therefore natural to wonder what occurs if these 
hypotheses are relaxed. 

The possibility of dimer dissociation introduces new time scales: the average lifetimes of all unstable j-clusters 
(2 < j < i* + 1). Within the framework presented in this paper, this problem is mapped into the random walk 
of a particle in a suitable high-dimensional space with, in general, a spatially varying diffusion coefficient. This 
inhomogeneity reflects the fact that unstable j-clusters diffuse and break up with rates different from the single 
adatom diffusion coefficient. In many cases the full solution is therefore beyond reach, even numerically (unless i* 
and L are very small) . 

However in the simplest cases our approach may still be fruitful. Let us consider for example i* = 2 and d = 1: 
three particles must meet (in the same lattice site) in order to nucleate a stable trimer. When two particles meet they 
form a dimer that dissociates after a typical time Tdis- If Tdis is much smaller than all other time scales, two adatoms 
diffuse as they were noninteracting, the 3d walker diffuses isotropically and we must just consider its irreversible 
passage along the diagonal (xi — xi — %■&). The same applies for generic d and i* as long as dissociation times of 
unstable clusters are small. This case is also of interest to test recent scaling approaches ]12] ] valid in the same limit 
(Tdis -> 0). 

The second natural extension of the present work consists in considering nucleation on top of fractal islands instead 
of compact ones. The framework of our method keeps unchanged Jl|. 

A further extension is to take into account the possibility of re-evaporation of deposited particles. 



APPENDIX A: ASYMPTOTIC BEHAVIORS OF THE TEMPORAL DISTRIBUTION IN D = 1. 

In this appendix we present some detailed results for the one-dimensional case with zero ES barriers and both 
adatoms having a uniform initial distribution. This is of course not physically sensible, since the effective initial 
distribution of the second adatom has a parabolic form for £ ES — 0. However, contrary to the physically sensible case, 
the evaluation of the coefficients B^- is not difficult and this allows an explicit analytic evaluation of the behavior of 
Q(t), which will differ from the realistic one only in the prefactors. 

A simple, although lengthy, evaluation of the coefficients leads to 
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We now want to calculate the temporal distribution 
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Q(t+1 >0) = 5^Pn-i,T.(t) 
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where the coefficients Ckj are given in Eq. (p9|). 

In the limit of large L, the coefficients Ckj are nonvanishing only for odd j — k (except for k — j, but Bkk = 0) and 
their value is 
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The dominant contribution to Q(£) comes from j = k — 1 or j = k + 1 
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We neglect the first term in both sums and approximate the sums according to the Euler-Maclaurin formula 
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fe=l 1/1 
Since in this case /l -C /1, we neglect /j,: 



rL_1 „ 2fc 2 + 2fc + l 

«fc — — — tttt- exp 



"9 6XP 



{-2k- l) 2 



5i / 7T 



t / 7T 



4 U + l 



(2fc 2 + 2k + 1) 



4 U + l 



L „ 2/c 2 - 2k + 1 
dfc (2fc-l) 2 GXP 



4 U + l 



(2fc 2 - 2fc + 1) 



(A15) 

(A16) 
(A17) 
(A18) 

(A19) 
(A20) 

(A21) 
(A22) 

If t <C (2L 2 /tt 2 ) each integral goes to and the first term prevails on the exponential; in the limit t ^> (2L 2 /tt 2 ) 



The integrals can be evaluated by considering for the integrand the limit for large k, and one obtains: 



Q(t) — 



irJt/2 



+ gjr^exp 



The upper limit of the integrals can be shifted to 00 



5t ( 7T 

4 \L+l 



5t ( it y 

4 \L+1 J 



the opposite is true. So, Q(t) ~ 8/(Lir 5 / 2 y/t/2) for t < (2L 2 /tt 2 ) and Q(t) ~ g^exp -f (x+l)' 
(2L 2 /tt 2 ). 

By summing Q(i) over t one obtains IT 



for t > 



IT 



<.2L 2 /tt 2 /■ 

~ / dtQ(t)= / 
Jo Jo 



2L 2 /iv 2 



dt- 
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L^l 2 sjtj2 7T 



r 7/2 



-0.58... 



(A23) 



f Email address: politi@fi.infn.it 
* Email address: castella@pil.phys.uniromal.it 
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